function [a b] = extract_curve_coeffs(e, de)
	% Little utility to calculate the coeffs in 'dE ~ a / (E + b)' given at
  % least two points on a dE-E curve, for cutting in dE-E plots.
	%
	%  [a b] = extract_curve_coeffs(e, de)
	%  eg. [a b] = extract_curve_coeffs([2000 12000], [4000 2000])

  % E = a / dE - b;
  % A = [1/dE -1];
  % y = E;
  % x = [a; b];

  A = [1./de' -ones(length(e), 1)];
  x = A \ e';
  a = x(1);
  b = x(2);
%	b = (de(1) / de(2) * e(1) - e(2)) / (1 - de(1) / de(2));
%	a = de(1) * (e(1) + b);
end
